Axion as a cold dark matter candidate: low-mass case 
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Axion as a coherently oscillating scalar field is known to behave as a cold dark matter in all 
cosmologically relevant scales. For conventional axion mass with 10 -5 eV, the axion reveals a 
characteristic damping behavior in the evolution of density perturbations on scales smaller than 
the solar system size. The damping scale is inversely proportional to the square-root of the axion 
mass. We show that the axion mass smaller than 10 -24 eV induces a significant damping in the 
baryonic density power spectrum in cosmologically relevant scales, thus deviating from the cold dark 
matter in the scale smaller than the axion Jeans scale. With such a small mass, however, our basic 
assumption about the coherently oscillating scalar field is broken in the early universe. This problem 
is shared by other dark matter models based on the Bose-Einstein condensate and the ultra-light 
scalar field. We introduce a simple model to avoid this problem by introducing evolving axion mass 
in the early universe, and present observational effects of present-day low-mass axion on the baryon 
density power spectrum, the cosmic microwave background radiation (CMB) temperature power 
spectrum, and the growth rate of baryon density perturbation. In our low-mass axion model we 
have a characteristic small-scale cutoff in the baryon density power spectrum below the axion Jeans 
scale. The small-scale deviations from the cold dark matter model in both matter and CMB power 
spectra clearly differ from the ones expected in the cold dark matter model mixed with the massive 
neutrinos as a hot dark matter component. 

PACS numbers: 98.80.-k, 95.35.+d 



I. INTRODUCTION 

Observations of type la supernovae, cosmic microwave 
background (CMB) radiation anisotropy, and large-scale 
structure of galaxies have revealed that our present uni- 
verse is dominated by two dark components [l|. Based 
on the Fricdmann world model, the energy content of the 
present universe is occupied by 73% of dark energy and 
23% of dark matter, whose nature is the fundamental 
mystery of the present day theoretical physics and cos- 
mology (see Ref . Q for recent cosmological interpretation 
of astronomical observations). 

The dark matter is often modeled as pressureless, non- 
relativistic particles, which is known as cold dark mat- 
ter (CDM). Although the CDM model has been greatly 
successful in explaining many cosmological observations, 
there remain certain conflicts at the galactic scales. First, 
the CDM-based simulations of galactic halos show den- 
sity profiles with a central cusp at kpc scale, while obser- 
vations indicate constant density cores [|| . Secondly, the 
CDM model predicts overpopulation of low mass halos 
within a galaxy or a galaxy group, which is an order of 
magnitude larger than the present observations [J. 

These CDM problems at galactic scales may be over- 
come by introducing a suppression of small-scale power 
in the density perturbation. There have been efforts to 
explain the CDM conflicts using the dark matter models 
with the Bose-Einstein condensate (BEC) and the ultra- 
light scalar field. 

In the BEC dark matter model the bosonic particles 
occupy the lowest quantum state of the external potential 



and can be described as a coherent matter wave of macro- 
scopic size due to Bose enhancement [oT-floj. Along this 
line there is fuzzy cold dark matter model proposed by 
Hu ct al. [TTJ] which is a case of free particles ignoring the 
self-interaction terms in the BEC based on the general- 
ized dark matter approach [l2| . In this model the coher- 
ent wave of ultra-light dark matter with m ~ 10 -22 eV 
can suppress the kpc-scale cusp in the dark matter ha- 
los and reduce the abundance of low mass halos. Recent 
investigations on the growth of perturbation in the uni- 
verse with BEC dark matter can be found in Refs. [9L [Tol| . 
The BEC dark matter has a close connection to the ax- 
ion dark matter model. Recently, Sikivic and Yang [l3j 
showed that the cold dark matter axions form a BEC 
through their self-interactions and gravitational interac- 
tions (see also Ref. [3). 

Another possibility is the ultra-light scalar field dark 
matter model [l5l - l24 | (see Ref. [2{| for a review). Based 
on a specific scalar field potential with appropriate initial 
conditions the evolution of scalar field mimics the prop- 
erty of the CDM fluid. For example, authors of Refs. 
[13 [3 proposed cosh-like potential V = Vq (cosh Ac/> — 1) 
as a good candidate of dark matter (see also Ref. [2fj|). 
Authors of Refs. [U, [22| investigated the structure for- 
mation in the scalar-field dark matter model with V ~ 
im 2 (/) 2 + iA0 4 by using the fluid and the field approaches. 
Authors of Refs. [HI, 13 studied the effect of ultra-light 
scalar field with V = ^m 2 <f> 2 on the growth of structure 
of the universe, in a situation where only a substantial 
fraction of ultra-light scalar field dark matter constitutes 
the total dark matter including CDM. 
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It is well known that the axion as a coherently os- 
cillating scalar field behaves as CDM in cosmologically 
relevant scales (26|. Khlopov et al. [27j investigated the 
growth of perturbations due to gravitational instability in 
a universe dominated by a scalar axion field. The quan- 
tum description of cosmological axion perturbations was 
given by Nambu and Sasaki (28|, while the full relativis- 
tic description based on the linear perturbation theory 
was presented in Refs. (29l - [3l| . For a typical mass with 
10~ 5 eV the behavior of the axion is equivalent to that 
of CDM on the cosmological scales including the super- 
horizon scale [3l| . 

In this work, we consider an axion dark matter model 
with extremely low mass (m < 10~ 22 cV) and investi- 
gate the effect of ultra-light axion on the baryon matter 
density and the CMB anisotropy power spectra and the 
perturbation growth. Our calculations are based on the 
full relativistic linear cosmological perturbation theory. 
Note that some of previous studies on BEC dark matter 
models relies on Newtonian cosmology and the relativis- 
tic description was given only in a limited fashion. Due 
to the common property of coherent oscillations, our de- 
scription for the low-mass axion is closely related with 
that from the BEC or ultra-light scalar field dark matter 
models. 

The structure of this paper is as follows. Section 
HT1 reviews the basic equations used to describe back- 
ground and perturbation evolutions based on the rel- 
ativistic linear perturbation theory, dealing with the 
multi-component fluids together with the minimally cou- 
pled scalar field. In Sec. [TIT] the basic properties of the 
axion dark matter are briefly described. A complete set 
of perturbation equations for the axion, baryon, and ra- 
diation fluids is presented, and initial conditions for the 
perturbation variables in the early radiation era are de- 
rived for two asymptotic limits. Our primary results arc 
summarized in Sec. IIV1 including baryon matter density 
and CMB power spectra and the growth of baryon den- 
sity perturbation for various values of extremely light ax- 
ion mass. We also compare the cosmological effects of 
light axion with the CDM model mixed with the mas- 
sive neutrinos as a hot dark matter component. We dis- 
cuss our results in Sec. [Vj Throughout this paper, we set 
c = 1 = h. 



II. BASIC EQUATIONS 

We consider scalar-type perturbations with the metric 



ds 2 = -(I + 2a)dt 2 - 2a/3 tCt dtdx° 



(l + 2tp)gj®+2 %a \p dx a dxP, (1) 



and introduce 



Here a(t) is the cosmic scale factor, a, f3, ip, 7 are metric 
perturbations, a dot and a comma indicate the time and 
spatial derivatives, respectively, and a vertical bar rep- 
resents the covariant derivative based on the comoving 

(3) 

three-space metric g a g . The energy-momentum tensor is 



^0° = ~M - Sfi, T^ = --(ji+p) v, a , 
T£ = (p + 6p)6%, 



(3) 



where we ignore the anisotropic stress; k is the wave num- 
ber with A = — k 2 in Fourier space. 

The background evolution is described by 



„, 8ttG K A 
H 2 = u + — . 

3 P a 2 3 

fi + 2,H{n+p) = 0, 



(4) 



where H = a/a is the Hubble parameter, K is the sign 
of spatial curvature, and A is the cosmological constant. 
In the multiple component case we have 



and 



j 3 



fii + 3H(^ +Pi) = 0, 



(5) 



(6) 



where we ignore direct interactions among fluids. For a 
minimally coupled scalar field we have the equation of 
motion 



<P + 3H ( f> + V 4> = 0, 
and the fluid quantities 



(7) 



M0 = ^ 2 + ^, P4> = \4?-V. (8) 
We need the following equations for the perturbations 

MM 



0, (9) 



5 fii + 3H(5ni + 8pt) - (fii +pi) [ K - 3Ha v t 

V a 
= 0, (10) 
I • k 

— [a 4 (/Zi +pi)vi] = - [(fii + p^ a + Spi] , (II) 



where 



X = a(f3 + aj), K = 3Ha-3(f> ^x- ( 2 ) 



j 3 

(fi+p)v = ^(Mj +Pj)Vj- 



(12) 
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In the presence of a scalar field we have the equation of 
motion 13311 



3HScj) 



2(pa + ( a + 6Ha + — X - 3<£ ] , (13) 



and the fluid quantities 

k ■ ■ ■ 

{fj-<j> +p<i>)v rjj = -<j)54>, 6(1$ = 4>5<j) - <j) 2 a + V^dcj), 



- 4> 2 a - V^8(j), 
with vanishing anisotropic stress. 



(14) 



where for the fluid quantities we have taken the time av- 
erage. Up to this point the perturbed quantities are spa- 
tially gauge invariant and the temporal gauge condition 
has not been taken yet [32j . 



A. Axion-comoving gauge 

As the temporal gauge (hypcrsurface or slicing) condi- 
tion we take the axion-comoving gauge [3lT | 



(Ma +Pa)v a = 0. 



Equation (|20|) gives 



(21) 



III. AXION 

We consider the axion in the presence of baryonic dust 
and radiation (photons and neutrinos). The axion is con- 
sidered as a minimally coupled scalar field with a poten- 
tial V = ^m 2 (j) 2 , where m is the axion mass. We strictly 
ignore H/m higher order term in our analysis. At the 
current epoch it is 



Ho 
m 



2.133 x 10" 



/ m 
UO" 5 



eV 



(15) 



where Hq = 100ft- kms _1 Mpc _1 is the Hubble expansion 
rate at present. 

In the axion model the equation of motion leads to the 
background solution as [US 



<j)(t) = a 3 / 2 [0+o sin(mt) + </>_ cos(mi)], 



(16) 



where 4>+o and 0-o ar e constant coefficients. We use the 
fluid formulation of axion. We have 

2 

Ma = £j(0+ O + 0-o), Pa = 0, (17) 

thus, the axion behaves as a zero-pressure fluid; for the 
fluid quantities the time-averaging has been applied for 
highly oscillating scalar field. For baryon, radiation (as 
a fluid), and axion, we have 



(J. = Mb + Mr + fj, a , p 



-U r . 



(18) 



To perturbed order, with an ansatz |29j 
6<f>(k, t) = 5<t>+{k, t) sin(mt) + 5<t>- (fe, t) cos(mt), (19) 
the fluid quantities in Eq. (|14p give [3l[ 
5[i a = a~ 3/2 m[ra(0 +o <50+ + 0_ o £0-), 



Sp a 

a 



-a 3/2 m(0 +o <50_ 



(fJ-a + Pa)v a = 



"3/2, 



b +0 6<t>- - 0- O <50+), (20) 



0+0 0-0 

and 

S Va 3/2 ^0+ Sp a 

= 2a ' a, = —a. 

Ma 0+0 Ma 

From Eqs. ([9]) and (fT0|) for the axion, we have 
k + 2Hk = A"KG{nb5b + 2fi r S r + fj, a S a ) 



5a = K, 



(22) 



(23) 





k 2 - 3K~ 


127rG ^ 6 + ^Mr) 


+ 2 



a, (24) 
(25) 



where Si = 6fii/[ii. For the baryon and radiation, under 
the fluid formulation, Eqs. (|TU)) and (fTTj) give 



8b = K — 3Ha Vb, 

a 

k 

v b + Hv b = -a, 
a 



SHa 



k ( ! r 
v r = - [ a + -8 r 

a V 4 



(26) 
(27) 
(28) 

(29) 



For a zero-pressure fluid (like a CDM) we have a = 
—Sp/n = 0. Now, in the case of axion, a can be de- 
termined in terms of 8 a using Eq. (|13p . From Eq. (|14| 
we have 



0+o 



a 3 / 2 k 2 



2m 2 a 2 0+o 
From these two relations, we have [3l[ 
k 2 1 

-8 a - 



Am 2 a 2 1 



k 2 

4m 2 a 2 



(30) 



(31) 



Equations (|2~4")) - (f2"9")) and (j3"Tj) constitutes the complete 
set for axion, baryon, and radiation. In a realistic situa- 
tion we have to use the Boltzmann equations for radiation 
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(photons and neutrinos) or th e tig ht-coupling approxima- 
tion for photon-baryon fluid [34 l35j. We note that the 
above equations are valid in the presence of K and A in 
the background. 

From Eqs. (JM]), (US]) and {gT]) we have 

5 a + 2H5 a - AttG (p a S a + ^bh + 2/i r 5 r ) 

k 2 - 3K 



12ttG ( n b + ~l-l r 



Am 2 a 2 1 



k 2 
4m 2 a 2 



0. 



(32) 



The equation is valid at all scales. Ignoring K , and for 
axion only, we obtain the axion Jeans scale, 



A./ = 



2na 



2n 



2-nGii a + y/^G 2 ^ 2 + 16TrG[i a m 2 



(33) 

In the limit of fc 2 /(4m 2 a 2 ) <C 1, with present Q a = 
8ir Gjj,„ /(3H 2 ) ~ 0.27, the current Jeans scale becomes 

[H M m M 



A.7 = 



G/i a0 m 5 



1/4 



50h 



-1/2 ^ 



10" 5 eV 



-1/2 



AU, 
(34) 



which is about the solar system size for m ~ 10 5 eV. 
For extremely low mass axion, we have 



X.j = 2Ah- 1/2 



10- 25 eV 



-1/2 



Mpc, 



(35) 



which is a cosmologically significant scale. On scales 
larger than the axion Jeans scale the axion behaves as 
the CDM. In the following we are interested in the ob- 
servable consequences of such a small mass axion as a 
dark matter candidate. 



B. Initial conditions 

In order to have initial conditions for perturbation vari- 
ables in the early radiation dominated era, we consider 
the baryon, radiation, and axion in a flat background. In 
the radiation dominated era (a cx t 1 / 2 ), Eqs. 
give 

■i 1 i 3 . / 3 fc 2 \ 

r 1 t ( 1 k2 \ c 2 i 2 . 3 

Sr + 2t Sr + [ ~¥ + 3^ ) 5r = ~3t 6a ~ I" + ^ 



(36) 



We consider asymptotic solutions in the large-scale limit, 
thus k/(aH) < 1. 



(I) For k 2 1 (m 2 a 2 ) <C 1, we have a = 0; this is the same 
as CDM. For 6 a oc t" we have (n-l)(2n-l)n(n+l) = 0, 
1, and the growing solution is 



thus n = 1, 2, 



<5a = h = -S r oc <. 



(37) 



This is the same as the adiabatic solution in the case of 
CDM. 

(II) For k 2 /(m 2 a 2 ) 3> 1, we have a = —S. For S a cx t n , 
we have (n — l)(n — 5)(n 2 +n+ |) =0, and the growing 
mode is 



5„ = —(5;, = — 5 r oc t. 
5 10 



(38) 



The condition k / (ma) 3> 1 together with k / (aH) <C 1 
demands H jm 3> 1, thus violating our basic assumption 
of the axion as a coherently oscillating scalar field. Al- 
though the numerical integration of the set of equations 
is possible, the equations are not valid in this case. In 
the small axion mass we are interested in, such a vio- 
lation is inevitable at some epoch in the early universe 
(Sec. HVA"| . To avoid this difficulty, as a simple solution, 
though ad hoc, we can introduce an evolving axion mass 
model (Sec. USUI . 



IV. OBSERVATIONAL EFFECTS 

In this section, we present our numerical calculations 
of background and perturbation evolutions in the axion 
dark matter model with extremely low mass. As a fidu- 
cial model, we use the flat ACDM model consistent with 
the Wilkinson Microwave Anisotropy Probe 7- year obser- 



vation (Qboh 2 



0.02260, fl c0 h 2 = 0.1123, fl A = 0.728, 



n s = 0.963, t = 0.087) including masslcss neutrinos with 
N v = 3.04 (see Table 14 of Ref. @j). As the dark matter 
component we consider the axion with Q ci q = Q c q. In 
order to calculate the matter and CMB power spectra, 
we solve a system composed of radiation, baryon, and 
axion in the axion-comoving gauge. The radiation com- 
ponents are handled by the Boltzmann equations while 
the photon-baryon fluid by the tight coupling approxi- 
mation. The detailed numerical methods are presented 
in Ref. M. 



A. Violation of axion-fiuid nature in the early era 

Due to the time dependence of H, our basic assump- 
tion of H/m 1 for axion as a coherently oscillating 
scalar field is violated in the early era especially for a 
low-mass axion. Top panels of Fig. [1] shows evolutions of 
H/m and k/(ma) for different values of axion mass and 
comoving scales. Despite such a violation of the basic 
axion-fluid assumption, here we present the baryon mat- 
ter density and CMB anisotropy power spectra for axion 
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FIG. 1: Top panels: evolution of H/m (left) and k/(ma) 
(right panel) as a function of scale factor a(t) for differ- 
ent choices of axion mass m and comoving wave number k. 
The left panel shows that for m = 10 _ eV, our basic as- 
sumption of H/m 1 is violated at a < 10~ 5,5 . On the 
right panel showing the evolution of k/(ma) for axion mass 
m = 10~ 24 eV, four different comoving scales from 10 -4 to 
0.1 /iMpc - have been chosen. Bottom panels: typical evolu- 
tion of density perturbations Si (i — a, b, c, 7) in the axion- 
comoving gauge. We considered the ACDM model (solid 
curves) and the A-axion-DM model with different choices of 
axion mass (dashed curves) . Here we consider two cases of ax- 
ion mass, m = 10 -24 eV (left) and 10 -22 eV (right panel). In 
both cases, compared with the cold dark matter in the ACDM 
model (S c ; blue solid curve), the axion density perturbation 
(5a', blue dashed curve) starts to evolve with lower amplitude; 
see Eq. p8[) . The behavior of 5 a slightly changes at a certain 
transition epoch which corresponds to when k/(ma) becomes 
unity. 



with a low mass. In section TlV Bl wc will introduce mod- 
els which avoid such a violation, and will show that the 
results are similar. 

Examples of evolution of individual density perturba- 
tion variables are shown in the bottom panels of Fig. Q] 
(see the caption for detailed descriptions), where the ini- 
tial conditions for k/(ma) 1 limit have been used for 
the axion mass and the comoving scales considered here. 
The violation of H/ m<l together with the large-scale 
assumption leads to the transition between k/(ma) <C 1 
and 3> 1 at the initial epoch (see Fig. [TJ top-right panel) . 
In our calculation we have imposed initial conditions for 
k/(ma) 3> 1 presented in Eq. ([55)1 because in most cases 
of comoving scale and axion mass considered the values of 
k/ (ma) are larger than unity. However, we have checked 
that the differences in power spectra obtained with dif- 
ferent sets of initial conditions in Eqs. ([3"7| and (|38l) arc 



A [h- 1 Mpc] 0=180*/; [degree] 
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FIG. 2: Baryonic matter (top-left) and CMB temperature, 
15-mode polarization, Ti?-cross power spectra (other panels) 
for axion mass 3 x 1CT 26 eV < m < 10~ 24 eV (colored curves ) . 
On the top-panels, the ratio of the power spectrum amplitude 
relative to the fiducial ACDM prediction (black curves) are 
shown to indicate the differences. As the mass decreases below 
10 -24 eV we have strong small-scale suppression of the mat- 
ter power spectrum due to the increase of axion Jeans scale. 
The CMB temperature power spectrum show small changes in 
higher multipoles corresponding to the scales below the axion 
Jeans scale. Compared with the baryon power spectrum, the 
change in the CMB power spectrum is minor. No noticeable 
changes appear in the polarization power spectra. 



negligibly small. We also note that if we apply the initial 
conditions depending on the value of k/(ma) it induces 
a significant discontinuity in the final power spectra. 

The matter and CMB power spectra expected in the 
low-mass axion model with a mass range 3 x 10 -26 eV < 
m < 10 -24 eV are shown in Fig. [51 As expected the 
small axion mass induces a significant damping in the 
baryon matter density power spectrum. As the axion 
mass decreases, the axion Jeans scale increases. From 
Eq. ([33)) we can see that for m ~ 1CP 25 eV we have 
the current Jeans scale Xj ~ 2.4 Mpc. However, the 
damping of the matter power spectrum appears at the 
comoving scale larger than the current Jeans scale. Note 
that the critical scale where the damping occurs is deter- 
mined by the Jeans scale at the radiation-matter equal- 
ity. Since the comoving Jeans wave number scales with 
the scale factor as kj cx a 1 / 4 during the matter era, 
we expect k Jcq m 0.15 /iMpc" 1 (Aj cq « 40 /j _1 Mpc) 
for Man = 0.27 and a cq ~ 3 x 10~ 4 (see Fig. HJ, see 
also [11[. Similarly, the CMB power spectra shows the 
deviation at scales smaller than the axion Jeans scale, 
thus at very high multipoles; for m = 3 x 10~ 26 eV, 
the damping of matter power spectrum occurs around 
a t kd > 0.1 hMpc^ 1 , which translates into the angu- 
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log a(t) 



FIG. 3: Evolution of H/m (left) and k/(ma) (right panels) 
as a function of scale factor a(t) for evolving axion mass model 
[Eq. (OH)]. Here o tr = 7xl(T 5 has been adopted and the axion 
mass ranges over 3 x 10 -26 eV < mo < 10 -5 eV. On the right 
panel the variation of k/(ma) is shown for a fixed comoving 
scale k = 0.1 /iMpc -1 and for axion mass 3 x 1CP 26 < mo < 
10 -20 eV. For comparison, the case of constant axion mass 
with m = 3 x 10 -26 eV is indicated by red dashed curves. 



lar scales i = r^cckd > 1000 (with rd oc ~ 10 4 /i _1 Mpc 
the distance to the decoupling surface) ; the deviation in- 
creases up to 15% at I ~ 3000 compared to the ACDM 
model (see Refs. [1, HI for predictions of the CMB 
anisotropy power spectra based on the BEC and the 
ultra-light scalar field dark matter models). 



B. Evolving axion mass model 

We have observed that the condition for axion as a 
rapidly oscillating scalar field is violated in the early uni- 
verse. As a simple solution to avoid this problem, we 
can introduce the evolving axion mass in the early uni- 
verse. In this way we can investigate the role of constant- 
low-mass axion in the later evolution of cosmic struc- 
tures within the axion-fluid formalism. One simple mass- 
evolution model is to assume that the axion mass varies 
in proportional to the Hubble parameter before a transi- 
tion epoch Otr and becomes constant after the transition. 
In the early epoch of radiation dominated era, the Hubble 
parameter varies with a~ 2 . Thus, assuming the transi- 
tion occurs in the radiation era, the evolving axion mass 
model can be written as 



m(a) 



m (a tr / a) 2 a < a tr , 
mo a > atr, 



(39) 



where mo is the constant axion mass after the tran- 
sition a tr . Alternatively, we may introduce a evolv- 
ing axion mass model with smooth variation with e.g., 
m(a) = mo[l + (an/a) 2 ]. Although not shown here, the 
results are quite similar to those from the simple evolving 
mass model adopted here. For the results shown in this 
section, we adopt a tr = 7 x 10~ 5 , which has been set from 
the condition that H/m « 0.1 at a < an for axion mass 
mo = 3 x 10 -26 eV, the lowest axion mass considered in 
this work. Figure [3] shows evolution of H/m and k/(ma) 



for different axion mass that varies with this functional 
form. 



C. Effects of small-scale damping on power spectra 
and growth rate 

Here we present the baryonic matter and the CMB 
anisotropy power spectra, evolution of perturbation vari- 
ables, and the growth factor based on the simple evolv- 
ing axion mass model in Eq. (|3"9"|) . The results for the 
power spectra are shown in Fig. U (solid curves) , where 
the power spectra for the constant axion mass (dotted 
curves; the same as in Fig. [2]) arc shown for comparison. 

Compared to the cases of constant axion mass, the 
power spectra for varying axion mass shows a small dif- 
ference: the damping of the power spectrum is slightly 
less significant as expected. The damping in the mat- 
ter power spectrum depends on the transition epoch an', 
more recent transition gives smaller damping. Compared 
with the ACDM fiducial model, the deviation in the CMB 
anisotropy power spectrum at high multipole I is now to 
the negative direction. Note that in the case of con- 
stant axion mass the deviation appears in the positive 
direction. In both cases, the magnitude of the maximum 
deviations are very similar to each other. 

In order to determine the low-mass limit of axion we 
present the baryon density power spectrum in Fig. [SJ 
which shows the substantial small-scale cut-off in the 
power spectrum for m < 10~ 22 eV. For m > 10~ 22 eV 
the matter power spectrum is almost identical to the 
CDM case on the scales considered (A < 3 /i _1 Mpc). 

Figure [6] shows the evolutions of density perturba- 
tion variables for photon, baryon, and dark matter in 
the evolving axion mass model at small comoving scales 
(k = 0.2 and 0.5 /iMpc -1 ) and for two extremely small 
axion mass (m = 3 x 10~ 26 and 2 x 10~ 25 eV). In all 
cases, the initial conditions for k/(ma) <C 1 limit have 
been used. 

In Fig. [7] we show the growth factor g = 5b /a (normal- 
ized to unity at present) at comoving scales k = 0.1, 
0.5 /iMpc -1 and for various axion mass values. At 
k = 0.5 /iMpc -1 , noticeable positive deviations from 
the ACDM prediction are seen. Despite the significant 
suppression of overall perturbation growth, the recent 
perturbation growth factor is larger in the low-mass ax- 
ion model, especially at the small comoving scales. This 
can be interpreted as the more growth in the axion mass 
model (see Refs. BGSIHLHi)- 



D. Mixed cold and hot dark matter model 

The warm dark matter model or the mixture of the 
cold and hot dark matter model are the alternative can- 
didates to the CDM model with small-scale suppression 
in the matter power spectra. It would be interesting to 
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FIG. 4: Baryonic matter (left) and CMB temperature anisotropy (right) power spectra for evolving axion mass model with 
tstr = 7 x 10 for different axion masses (solid curves). Dotted curves represent the similar power spectra obtained when the 
constant axion mass is assumed which are the same as those in top panels of Fig. [2] The bottom panels show power ratios 
relative to the fiducial ACDM model prediction. In the CMB power spectrum the deviations at high multipoles (£ > 1000) 
show opposite trends depending on the constant mass and the evolving mass models. 
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FIG. 5: Baryonic matter density power spectra expected 
in the evolving axion mass model with 10 -25 eV < mo < 
10 -22 eV (solid curves). Now the range of comoving wave 
number has been extended to the higher values (smaller 
scales) around k — 2 ftMpc -1 . Similarly to Fig. 01 the 
power spectra for the constant axion mass are shown as dot- 
ted curves. The bottom panel shows the ratio of power rela- 
tive to the ACDM prediction, indicating the power spectrum 
damping on small scales. 




-10 -8 -6 -4 -2 -10 -8 -6 -4 -2 

log a(t) log a(t) 

FIG. 6: Evolution of density perturbation variables for pho- 
ton, baryon, dark matter for evolving axion mass model, Eq. 
(|39[) with a tr = 7x 10~ 5 . Shown are typical evolutions of den- 
sity perturbations Si (i = a, b, c, 7) in the axion-comoving 
gauge for m = 2 x 10 -25 , 3 x 10 -26 eV and k = 0.2, 
0.5 /iMpc - (dashed curves). We also considered the ACDM 
model whose density perturbations are shown as solid curves 
in each panel. Compare the results with those in Fig. [1] (bot- 
tom panels) 
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FIG. 7: Evolution of the normalized perturbation growth fac- 
tor g = (St/a) at comoving scales of k = 0.1, 0.5 /iMpc -1 
for different axion masses (colored curves). Here we use an 
evolving- axion- mass model with a tT — 7 x 10~° . Black curves 
represent the growth factor of the fiducial ACDM model. 



X [h" 1 Mpc] 0=180"// [degree] 
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FIG. 8: Matter and CMB temperature anisotropy power spec- 
tra expected in a model containing small fraction of massive 
neutrino component, where we have used the CAMB soft- 
ware [37}|. Here we consider massive neutrinos (N v = 3.04) 
whose mass ranges from m v = 0.154 eV (fi„o = 0.01; red) 
to 0.769 eV (fi„o = 0.05; violet curves), with a relation 
(Q. v o + fl c o)h 2 = 0.1123. Black curves represent the power 
spectrum of the fiducial ACDM model with massless neutri- 
nos. The curves in the bottom panels indicate the ratios of 
powers relative to the ACDM model prediction. 

compare the small-scale suppression of our low-mass ax- 
ion model with those models. In Fig. [5] we present the 
case of massive neutrinos contributing as the hot dark 
matter added to the dominant CDM model. The results 
are similar to those found in Ref. [3t| . 

Although the low-mass axion dark matter model shows 
a sharp damping in the baryonic matter power spectrum, 
the CDM model with a small fraction of massive neutri- 



nos gives a damping of power with almost constant fac- 
tor at k > 0.1 /iMpc -1 . Besides, the CMB anisotropy 
power spectrum is very sensitive to the massive neutrino 
contribution; the massive neutrinos affect all the acous- 
tic oscillatory features at intermediate and high angular 
scales (£ > 100). On the other hand, the behavior of the 
axion dark matter model is quite different, and is almost 
insensitive to the axion mass except at the higher multi- 
poles. In the low-mass axion case the changes from the 
CDM in both power spectra occur only at scales smaller 
than the axion Jeans scale. 



V. DISCUSSION 

In this work, we have studied effects of extremely low- 
mass (to < 10 -24 eV) axion on the baryon matter density 
and the CMB anisotropy power spectra, and the pertur- 
bation growth based on the full relativistic linear per- 
turbation analysis. With a low mass, however, the basic 
assumption about the coherently oscillating scalar field 
is inevitably broken in the early universe (H/m 3> 1). 
We have introduced the simple evolving-axion-mass-in- 
the-early-universe model to avoid this problem. 

We showed that axion mass smaller than 10~ 24 eV in- 
duces the characteristic significant damping in the baryon 
density power spectrum on scales smaller than the axion 
Jeans scale, and changes in the higher multipole in the 
CMB anisotropy power spectra. Except for small changes 
in the higher multipoles (£ > 1000) corresponding to the 
scales smaller than the axion Jeans scale, the CMB power 
spectrum remains the same as the CDM case. The CDM 
nature is also preserved in the baryon matter power spec- 
trum above the axion Jeans scale. We showed that the 
small-scale damping nature of our low-mass axion model 
differs from the one expected in the CDM model mixed 
with the massive neutrinos as a hot dark matter compo- 
nent. 

Whether the small-scale damping of the baryon mat- 
ter density power spectrum can help alleviating the ex- 
cess small-scale clustering problem of the CDM model 
requires further studies in the nonlinear clustering prop- 
erties of the light mass axion. 
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